Computing Aα, log(A) and Related Matrix Functions by Contour Integrals
نویسندگان
چکیده
New methods are proposed for the numerical evaluation of f(A) or f(A)b, where f(A) is a function such as A1/2 or log(A) with singularities in (−∞, 0] and A is a matrix with eigenvalues on or near (0,∞). The methods are based on combining contour integrals evaluated by the periodic trapezoid rule with conformal maps involving Jacobi elliptic functions. The convergence is geometric, so that the computation of f(A)b is typically reduced to one or two dozen linear system solves, which can be carried out in parallel.
منابع مشابه
Computing AAlpha, log(A), and Related Matrix Functions by Contour Integrals
New methods are proposed for the numerical evaluation of f (A) or f (A)b, where f (A) is a function such as A 1/2 or log(A) with singularities in (−∞, 0] and A is a matrix with eigenvalues on or near (0, ∞). The methods are based on combining contour integrals evaluated by the periodic trapezoid rule with conformal maps involving Jacobi elliptic functions. The convergence is geometric, so that ...
متن کاملEvaluating Matrix Functions for Exponential Integrators via Carathéodory-fejér Approximation and Contour Integrals
Abstract. Among the fastest methods for solving stiff PDE are exponential integrators, which require the evaluation of , where is a negative semidefinite matrix and is the exponential function or one of the related “ functions” such as . Building on previous work by Trefethen and Gutknecht, Minchev, and Lu, we propose two methods for the fast evaluation of that are especially useful when shifte...
متن کاملThree-Dimensional Interfacial Green’s Function for Exponentially Graded Transversely Isotropic Bi-Materials
By virtue of a complete set of two displacement potentials, an analytical derivation of the elastostatic Green’s functions of an exponentially graded transversely isotropic bi-material full-space was presented. Three-dimensional point-load Green’s functions for stresses and displacements were given in line-integral representations. The formulation included a complete set of transformed stress-p...
متن کاملComputing the Matrix Geometric Mean of Two HPD Matrices: A Stable Iterative Method
A new iteration scheme for computing the sign of a matrix which has no pure imaginary eigenvalues is presented. Then, by applying a well-known identity in matrix functions theory, an algorithm for computing the geometric mean of two Hermitian positive definite matrices is constructed. Moreover, another efficient algorithm for this purpose is derived free from the computation of principal matrix...
متن کاملUnderground contour (UGC) mapping using potential field, well log and comparing with seismic interpretation in Lavarestan area
Coastal Fars gravimetry project in Fars province was carried out to find the buried salt domes and to determine characteristics of faults in this area. The Lavarestan structure was covered by 4203 gravimetry stations in a regular grid of 1000*250 m. Depth structural model of this anticline made in previous studies was based on geological evidences and structural geology measurements. In order t...
متن کامل